
00501 6044 


NASA Contractor Report 165793 




NASA-CR- 165793 

ocoHoC E' 


BIAXIAL TESTS OF FLAT GRAPHITE/EPOXY 
LAMINATES 


Harold Liebowitz and Douglas L. Jones 


THE GEORGE WASHINGTON UNIVERSITY 
School of Engineering and Applied Science 
Washington, DC 20052 


Grant NSG-1289 
October 1981 


bin i 


NOV 1 8 1981 


L'".r.’GLEY 

l:g 


^ A R F ♦ r,r r,i j p 
iRAfi’Y, NASA 

"CM. v!rg:n:a 


IVV\SA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



SUMMARY 


This report summarizes the research performed under NASA 
Grant NSG1289 Improved Fracture Toughness Determination for 
Aircraft Structural Materials, over the life of the Grant, : 
February 1, 1976, to December 31, 1980. The Technical Officers 
of this Grant were Drs. Wolf Fiber and G. L. Roderick of the 
NASA Langley Research Center. 

The primary purpose of this research program was to ex- 
amine analytically and experimentally the influence of biaxially 
applied loads on the strength of composite materials containing 
holes. The analysis was performed through the development of a 
three-dimensional, finite-element computer program that is cap- 
able of considering three types of damage: fiber breakage, 

delamination, and matrix failure. Realistic failure criteria 
were established for each of the failure modes and the influence 
of biaxial loading on damage accumulation under monotonically 
increasing loading was examined in detail. 

The experimental portion of this program encompassed both 
static and fatigue testing of specially designed biaxial speci- 
mens containing central holes. Initially, static tests were 
performed to obtain an understanding of the influence of bi- 
axial loads on the fracture strength of composite materials and 
to provide correlation with the analytical predictions. The 
predicted distributions and types of damage were in reasonable 
agreement with the experimental results. Then a number of 
fatigue tests were performed to determine the influence of 
cyclic biaxial loads on the fatigue life and residual strength 
of several composite laminates. 
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LIST OF SYMBOLS 
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(|)^ = the angle between X-axis and the fiber orientation in 
the i-th layer 

a., = stress tensor 


= strain tensor 

^iimn ~ niaterial constants for linear elastic anisotropic 
solid 

(L,T,z) = coordinate system; L-axis is parallel to the fiber 
orientation; z-axis is perpendicular to the 
laminate 

= Young's moduli 


^LT’^Tz ~ Poisson's ratios 
^LT’*^Tz ~ shear moduli 


d. . 
ij 


material constants for 
expressed in (L,T,z) 


elastic solid with 
coordinate system 


axis symmetry 


d. . 
ij 


material constants for 
expressed in (x,y,z) 


elastic solid with 
coordinate system 


axis symmetry 


~ displacements in x,y,z direction, respectively 


t^ = stress vector 


= unit normal vector of a surface 

(r,0,z) = cylindrical coordinate system 
Tq = radius of the centered hole 

[6] = nodal point displacement vector 
= shape functions 

- nondimensional coordinates for a three dimensional prism 
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[B] = matrix linking strains and nodal point displacement vector 

[J] = Jacobian matrix 

[F] = nodal point force vector 

[K] = stiffness matrix 

,o^ ,o^ ,Oj^^ ,Oj^^ , 0 ^^ = stress components in (L,T,z) coordinate 

system 

*^FN ’*^FS ’^MN *^MS ’^DN ’*^0 5 ^ failure strengths of a unidirectional 

fiber-reinforced composite layer 

R^,R 2 »...Rg = ratios between stresses and failure strengths 


R[I] = damage ratio of the I-th element 

R„,^ = maximum damage ratio of the entire laminate 

m 3. X 


0 = applied stress 
k = biaxial load factor 

2L,2W,2h = Total length, width, height of the laminate 



INTRODUCTION 


The research program. Improved Fracture Toughness Deter- 
mination for Aircraft Structural Materials, was initiated for 
the purpose of studing several problems of interest to NASA. 

The grant, NSG1289, was sponsored through the NASA Langley Re- 
search Center (LaRC) over the period February 1, 1976, to 
December 31, 1980. 

The initial task of this grant was to perform a detailed 
assessment of several nonlinear fracture toughness parameters, 
with the emphasis to be placed on the two-parameter Newman ’ 
method. This assessment was completed and the results reported 
to LaRC in September, 1976. 

The second and most important task was to evaluate both 
analytically and experimentally the effect of biaxial applied 
loads on composite laminates containing holes. In order to 
perform this research it was necessary to develop a compre- 
hensive three-dimensional, finite-element program and an exr 
perimental facility capable of applying substantial biaxial 
loads on composite specimens. Both of these capabilities were 
established and subsequently employed in this research program. 
The primary emphasis of this program was placed on multilayer 
composite materials in which the individual plies were oriented 
at 0,±45 and 90 degrees from an axis of symmetry. This con- 
figuration is approximately isotropic in its material behavior 
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and because o£ this property has been widely employed in 
aerospace structures. In addition, some testing has been 
performed on [0,±45,0]j^g laminates (n = 1,2). 

Both the analytical and experimental tasks emphasized 
specimens containing central holes, since they provide free 
surfaces at which the damage is expected to initiate. The 
diameter of the central hole was not varied in the analysis, 
although the hole diameter-to-specimen width ratio (2rQ/2w) 
could easily be varied in subsequent computations. Since 
the test specimen size was fixed, some experiments were 
performed with different hole diameters so that the influence 
of hole size could also be examined. 

The emphasis in the finite-element analysis was placed 
on monotonically increasing loading, because of the difficulty 
and expense of applying three-dimensional, finite-element 
programs to fatigue problems. Therefore, the initial experi- 
ments were also performed statically for the purpose of ob- 
taining comparisons with the experimental results. In addi- 
tion, a substantial number of fatigue tests were performed, 
so that an understanding could be obtained of biaxial load 
effects on the fatigue life and residual strength of composite 
specimens with holes. The principal results of the analytical 
and experimental results are discussed in greater detail in 
subsequent sections of this report. 
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FINITE-ELEMENT ANALYSIS OF LAYERED COMPOSITE MATERIALS 


One of the two primary objectives of this grant was the 
development of a three-dimensional, finite -element computer 
program to analyze damage accumulation in composite laminates. 
This program has been developed and is capable of: 

(i) providing a detailed stress distribution in a 

structural member made of a layered fiber-reinforced 
composite material containing a round hole and sub- 
jected to biaxial applied loading; 

(ii) identifying the shape and extent of the damage zone 
and the mode of failure; and 
iiii) analyzing the problem of accumulated damage and 

progressive failure for any given failure criterion. 

The principal characteristics of this program are outlined in 
the subsections below. 

Problem Description 

The specimen geometry is that of a rectangular plate 
made from a layered composite material containing a centered 
circular hole and subjected to in-plane biaxial loading. 

Each layer is a unidirectional fiber reinforced composite 
material oriented at some angle (j> with respect to an axis 
of symmetry (the x-axis) . The input data for a graphite- 
epoxy laminate are given below. 
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1. Material Constants 

E, = 134.58 X 10^ MPa 
L 

= 9.68 X 10^ MPa 
= 0.264 
v^2 = 0.435 

= 4.48 X 10^ MPa 

Li 1 

2. Ply Configuration 

[0,45,-45,90]^ 

3. Specimen Parameters 

2L = 152.4 mm 
2W = 152.4 mm 
2h = 1.07 mm 
r^ = 12.7 mm 

4. Finite-Element Mesh 

Number of nodal points = 1750 
Number of degrees of freedom = 5250 
Number of elements = 1248 

The coordinate directions are defined in Fig. 1. 

Linear Elastic Theory 


Each layer is modeled as a linear elastic continuum 
having an axis of symmetry about its fiber orientation. 
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Until de lamination occurs the layers are assumed to be 
perfectly bonded together. No adhesive layer is included 
in this study. 

Finite-Element Method 

The finite-element mesh consists of a family of three- 
dimensional prisms, each having eight nodes. Each node is 
assumed to have three degrees of freedom of motion. The 
symmetry and periodicity of this geometry have been taken 
into consideration. The initial results of this analysis have 
been published by Lee [1] . 

Failure Criteria 

It was assumed that only three types of damage can occur 
in this laminate: fiber breakage, debonding, and matrix 

failure. Failure criteria (strengths) associated with each 
mode of failure are given below 

= 1241 MPa 

= 82.7 MPa 

o„„ =55.2 MPa 
MN 

o„_ =82.7 MPa 
MS 

o^.. = 55.2 MPa 
DN 

( 

Opg = 89.6 MPa 


The stress components are defined in Fig. 1. 
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The Numerical Procedure 

An iterative procedure has been devised to identify 
the damaged zone and the mode of failure. After this has 
been done, the appropriate constraining forces are reduced 
to zero and the local stiffness matrix of the damaged zone 
is modified according to the mode of failure. Then the stress 
distribution is recalculated to determine whether the damage 
is progressing. When the damage is determined to be 
propagating the ultimate strength of the laminate is then 
assumed to be known as a function of the specific biaxial 
load ratio. 

The detailed formulation, analysis, and the numerical 
results are contained in a paper entitled, "Three Dimensional 
Finite-Element Analysis of Damage Accumulation in Composite 
Laminate," which has been accepted for publication in 
Computers and Structures [2] . A copy of this paper is 
included as Appendix A of this report. The principal 
conclusions of this paper are listed below. 

Cl) The initial stress distribution of near the edge 
of the hole determines the location of the initial 
damage zone and its mode of failure as well as the 
applied stress necessary to cause initial damage. 

(2) The initial stress distribution of near the 

edge of the hole predicts the final path of fracture. 

(3) No delamination was identified at any time in the 
failure process. 
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EXPERIMENTAL PROGRAM 


The primary objective o£ the experimental portion of 
this research program has been to establish a data base 
ivhich vrould serve to identify the effects of biaxially 
applied loads on the failure characteristics of composite 
materials. The experimental research program required the 
following initial steps to get the test system into operation 

(i) the design and fabrication of a biaxial testing 
system, 

(ii) the design of a suitable test specimen or specimens, 
and 

(iii) the performance of some initial tests on metal 
specimens to check out the test system. 

These tasks were completed early in the research program as 
indicated in the early progress reports. The test system 
developed for these tests is illustrated in Fig. 2 and the 
general specimen geometry is shown in Fig. 3. A more 
detailed description of the system development was provided 
in an earlier progress report [3]. In addition, the results 
of a number of tests on metals and polymers have been 
reported [4] . 

After the test system was determined to be functioning 
properly, five test series were conducted in the following 


sequence . 
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(1) Four initial static strength tests were performed 
on graphite/epoxy specimens, having [0,±45,90]^ 
and [0,±45,90]2g configurations, in which central 
hole diameters of 12.70 and 50.80 mm were 
utilized. The primary purpose of these tests was 
to determine whether variations in hole size in- 
fluenced either the strength or the magnitude of 
the biaxial effect, for k = 0 and 1.0. 

(2) Seven additional static strength tests were per- 
formed on graphite/epoxy specimens having 
[0,±45,90] , [0,±45,90] and [0,±45,0]_ configura- 
tions as a further evaluation of biaxial and hole 
size effects. Center hole diameters of 25.4 and 
50.8 mm and biaxialities of k = 0, 1/2, and 1.0 
were used in these tests. 

(3) After evaluation of the breaking strengths of the 

first eleven specimens, a series of seven fatigue 
tests were conducted on [0,±45,0] , and 
[0,±45,0]2g laminates. The test conditions for 
these specimens were: center hole diameters - 

25.4 and 50.8 mm; loading frequency - lOHz in- 
phase cyclic loading; and biaxiality ratio - 

k = 0, 0.25 and 0.50. 
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(4) A series o£ one static and five fatigue tests was 
performed on six additional specimens having a 
[0,±45,90]g configuration. The static test was 
performed for the purpose of comparing with 
earlier specimens having the same layup pattern. 

The test parameters for the fatigue tests were 
similar to the previous fatigue tests except that 
the biaxiality ratios were k = 0, 0.5, and 1.0. 

(5) A series of one static and five fatigue tests was 
performed on six specimens having the [0,±45,0]^ 
layup pattern. The initial static test in this 
series was different from the previous ones in 
that the loading rate was the same as that for the 
cyclic tests (lOHz) whereas the previous static 
tests had been performed slowly (for data-taking 
purposes) . The remaining fatigue tests were per- 
formed as in the earlier tests for this layup 
except that 50.8 mm diameter center holes were 
used for all specimens. 

After all of the test results had been evaluated, their 
variations with biaxiality ratio, hole diameter, thickness, 
etc. , were collected and examined. These results are included 
in this report, along with the conclusions drawn from them. 



10 


Test System Development 

The biaxial test system illustrated in Fig. 2 was de- 
veloped by fabricating a special, light-weight test frame of 
222,400 N (50,000 lb.) capacity and interfacing it with an 
existing 444,800 N (100,000 lb.) capacity test system used 
as the vertical axis. Both systems were of the servo- 
hydraulic type so that both static and dynamic loads could 
be applied to either axis. An additional control system was 
added that has the capability of controlling both axis and 
special equipment was added to eliminate interaction or 
’’cross- talk" between the two channels. Since only two 
hydraulic actuators were used in this system, the center 
of the test specimen would not remain stationary but would 
move in response to each actuator. In order to prevent 
this motion from introducing side loads onto the specimen 
the horizontal axis was suspended in its working position 
by elastic ropes. The deformation and pendulum action 
permitted by the ropes permitted slight movement in both 
the vertical and horizontal directions with the center of 
the specimen. Although the motion of the frame under rapid 
cycling transmits dynamic inertial loads to the specimen, 
the magnitudes of these loads were calculated for expected 
worst-case conditions and were found to be small enough to 
be neglected. 

A biaxial test specimen was designed as a modification 
of an existing specimen configuration used elsewhere [5-7]. A 
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photoelastic study was initiated for the purpose of maxi- 
mizing the size of the central uniformly-stressed work area. 

It was also necessary to design and purchase (or fabricate) 
grips to interface between the test system and the specimen. 
The specially designed grips that were purchased performed 
satisfactorily throughout the testing program. 

After assembly of the horizontal test system, it was 
operated for a brief check-out period and then interconnected 
with the existing vertical test system. The combined system 
was then checked out and found to operate as designed. This 
system design permitted phased operation of the two channels 
as well as completely independent operation. 

Experimental Procedure 

While the test system was being developed, the character- 
istics (material, layup, etc.) of the various test specimens 
were determined through discussions between the technical 
officer and one of the principal investigators. A typical 
specimen is shown in Fig. 4, where it is seen that a rein- 
forcing aluminum sheet was bonded to each side of the speci- 
men in the area around the four loading tabs. The specimens 
were manufactured by a separate vendor to the specifications 
agreed upon by the technical officer and principal investi- 
gator. 

While the specimens were being prepared, the test con- 
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ditions (hole diameter, biaxiality ratio, loading conditions, 
etc.) were also established for the first test series of 
four static tests. A total of thirty specimens, all made of 
5208-T300 graphite/epoxy lamina, were prepared and tested 
under this program, the characteristics and test conditions 
for these specimens are shown in Table 1. Biaxiality ratios 
were limited to the range 0 < k < 1.0 for two reasons: (i) 

the quasi-isotropic [0,±45,90]^^ configuration is symmetric 
with respect to the horizontal and vertical axes, so that the 
reciprocal biaxiality ratios, °° < k < 1.0 are also included 
in these tests, and (ii) the [0,±45,0]^^ configuration was 
limited to ratios of 0 < k < 0.5, because of its reduced 
strength in the 90° orientation. 

The specimen design shown in Fig. 4 created one type 
difficulty that had considerable influence on some of the 
test parameters. The problem was caused by the tendency of 
some specimens to fail by tearing off one or more of the 
loading tabs, resulting in the lack of failure through the 
test section. These tests could not provide valid results 
so certain modifications were introduced in'*‘o the specimen 
design to reduce the tendency for tab failures. The problem 
was particularly bothersome in the [0,±45,0]^^ specimens be- 
cause the tabs in the 90° orientation would fail under quite 
small biaxiality ratios. Because of this problem most of 
the later tests employed hole diameters of 50.8 mm. This 
modification was still inadequate for several tests so 
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additional modifications were introduced. The most success- 
ful modification was to machine away the reinforcing 
aluminum beyond where the individual tabs were separated. 

A specimen modified by these procedures is shown in Fig. 5 
and such modifications were employed in many of the later 
tests. The problem of tab tear-off was also greater for the 
fatigue tests than for the static tests. 

Test procedures unique to each test series are now dis- 
cussed in some detail so that the significance of the results 
may be more fully appreciated. 

Test Series (1) 

This test series was comprised of four static strength 
tests on the quasi-isotropic laminates involving variations 
in thickness (8 and 16 ply laminates), hole diameter, 

(D = 12.7 and 50.8 mm) and biaxiality ratio, k = 0 and 1.0. 

The loads were applied as in-phase ramp functions having 
300 second duration. The ramp functions were interrupted 
periodically to record data manually from digital multimeters. 
Five channels of data were recorded, including load and 
stroke on each axis and the hole opening displacement in the 
primary (largest load) direction. Load displacement curves 
in the primary direction were also recorded using an x-y 
recorder. All specimens failed through the center hole. 
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Test Series (2) 

Upon satisfactory completion of the first four tests, 
seven additional static tests were performed on both the 
quasi-isotropic and the [0,±45,0]2g configurations. The 
first two tests were performed to examine the influence of 
fabrication damage, which existed in a number of the speci- 
mens, around the central hole. The most extensively 
damaged specimen was 9-2 and the extent of the damage is 
indicated in Fig. 6a. Identical test conditions were applied 
to specimens 9-2 and 9-3 (which had negligible damage around 
the hole) , with the test results indicating that the initial 
damage did not influence the strength. The breaking strength 
was slightly greater for 9-2 and the final fracture occurred 
away from the initial damage, as seen in Fig. 6b. As a 
result of these tests, the existing initial damage was not 
taken into consideration in the subsequent testing. 

Specimens 5-3, 5-4 and 5-5 represented a test series in 
which only the biaxial load ratio was varied. The material 
was quasi- isotropic and biaxial ratios of 0, 0.5 and 1.0 
were selected. The tests were performed using in-phase 
ramp functions on each axis and all failures traversed 
through the center holes. 

Specimens 11-1 and 11-2 also represented a biaxial static 
strength series for the [0,±45,0]2g configuration. Biaxial 
ratios of k = 0 and 0.5 were employed for these tests and the 
failures traversed through the center holes. 
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Test Series (3) 

This test series was comprised o£ seven fatigue tests 
on the [0,±45,0], and [0,±45,0] configurations. The 
purpose of these tests was to determine whether damage 
accumulation due to cyclic loading at low biaxiality ratios 
would have any significant influence on the strength in the 
0® fiber orientation. Based on the static strength tests 
already performed on the [0,±45,0] configuration (Speci- 
mens 11-1 and 11-2) , the first specimen (11-3) in this 
series was fatigue cycled in uniaxial tension at 75 percent 
of the ultimate strength obtained from Specimen 11-1. The 
specimen had not failed after 10^ cycles so the residual 
strength was measured following the same procedures as for 
the static strength tests. The next specimen was fatigue 
cycled under the same load in the 0° direction but with 
k = 0.5 rather than 0. This specimen failed quickly by 
tearing off one of the horizontal tabs (90° direction) and 
thus did not represent a strictly valid biaxial comparison. 
The third specimen was cycled at the same primary load level 
as before but with k = 0.25. This specimen failed after 
approximately 440,000 cycles by pulling off one of the 
primary-axis tabs so it also could not represent a valid 
biaxial comparison. 

Because of the difficulties encountered in obtaining 
valid biaxial comparisons in the three previous tests, the 
next four tests (Specimens 7-1, 7-3, 7-4, and 7-5 in the 
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[0,±45,0] configuration) were employed primarily for deter- 
mining a modified specimen design capable of eliminating the 
tab failure problem. All four tests were performed with a 
biaxiality ratio of 0.25. The first specimen was tested 
without modification as a check on the necessity of speci- 
men modification for the eight-ply configuration. The 
applied cyclic loading for 7-1 was selected to be 80 per- 
cent of the strength of the corresponding sixteen-ply speci- 
mens and the biaxiality ratio was k = 0.25. This specimen 
sustained 10^ cycles after which the residual strength was 
evaluated and found to be significantly greater than the 
estimated original strength. This specimen failed through 
the center hole and along one edge so it was considered to 
represent a valid failure. The next specimen (7-3) was 
modified by increasing the center hole diameter from, 25. 4 to 
50.8 mm and was tested at a considerably higher cyclic load. 

It failed after approximately 3 X 10^ cycles by tearing off 
two adjacent loading tabs and so did not provide valid com- 
parisons. The next specimen (7-4) was modified by increasing 
the depth of the twenty slots between the individual loading 
tabs to the edge of the aluminum reinforcing plate. The slots 
were extended for two purposes: (i) to further reduced the net 

section of the specimen through the center hole, and (ii) to 
prevent the initial cracking of the aluminum from contri- 
buting to stress concentration effects. This specimen failed 
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during cyclic loading by tearing off two adjacent loading 
tabs. Some of the damage extended to the center hole but 
the failure did not pass through it. 

The final specimen in this series was also modified by 
increasing the central hole diameter to 50.8 mm and ex- 
tending the slot depth as was done for the previous test. 

This specimen also failed by pulling off two perpendicular 
tangs, although the 0° direction tab failure extended through 
the center hole. This test was thus considered to be mar- 
ginally successful in performing the desired objective. A 
photograph of the fractured specimen is shown in Fig. 7. 

Test Series (4) 

This test series was comprised of one static and five 
fatigue tests on specimens having the quasi-isotropic 
[0,±45,90]g configuration. The first test was a static 
test for the purpose of comparing the strength with prior 
tests of the same configuration. The test conditions for 
this specimen, 23-1, were identical to specimen 5-3 and the 
ultimate loads were in very good agreement with each other. 
Since there was no evidence that the quasi-isotropic speci- 
mens had a tendency for failure at the tabs, a fatigue test 
was performed at 80 percent of the ultimate (k = 0 in both 
cases) and the specimen failed at approximately 218,000 
cycles by tearing off one of the tabs in the primary 
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direction. It is possible that the failure was due to a 
momentary overload caused by a power failure in the labora- 
tory so the absolute value of the cyclic life was not 
reliable. Since the failure did not occur through the 
center hole, the test would not have been suitable for 
comparison purposes. For the remaining four specimens the 
center hole diameter was increased to 50.8 mm to increase 
the probability of center failures. The next three tests 
(23-3, 23-4, 23-6) represent the influence of biaxial 
loads (k = 0, 0.5 and 1.0) on the fatigue life and all 
three provided valid data. The final test, 23-5, was 
listed at a value of k = 0.5 and sustained 10^ cycles. The 
residual strength was tested and found to be quite high. 
However, this specimen failed by breaking off of two 
perpendicular tabs and the corner of the test area between 
them. Since the fracture did not go through the center 
hole the results were not considered to be valid. 

Test Series (5) 

This test series was comprised of one static strength 
test and five fatigue tests on specimens having the 
[0,±45,0]^ configuration. By the time that this test series 
was conducted, a PDP-11-03 computer had been added to the 
system for the purpose of rapid data collection. Therefore, 
it was decided that the static test would be performed at a 
rate of lOHz, rather than at the low rates previously 
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employed. It was hoped that a test of this type would pro- 
vide a more accurate measure of the ultimate strength and 
thus improve the accuracy of the percent of ultimate stress 
used in performing the fatigue tests. Also, because of the 
prior problems encountered with tab failure, two modifica- 
tions of these specimens were made: (t) the hole diameter 

was increased from 25.4 to 50.8 mm, and the aluminum 

reinforcing plate was machined away to approximately 6.3 mm 
beyond the slot tips without disturbing the composite 
material. These modifications were successful, since all 
six specimens failed through the center hole. A photograph 
of one of the modified specimens is shown in Fig. 5. Speci- 
mens 24-2 and 24-3 were loaded at 80 percent and 90 percent, 
respectively, of the ultimate strength determined from speci- 
men 24-1. Both of these specimens sustained 10^ cycles and 
provided high residual strengths. The biaxial ratio was 
0.25 for all three tests. 

The final three tests represented a series of three 
fatigue life tests in which the biaxial ratio had the values, 
k = 0, 0.1, and 0.25. All three specimens were loaded at 
least 90 percent of the static strength obtained from speci- 
men 24-1 and all failed in less than 10^ cycles. Thus, the 
results were not as conclusive as they would have been if 
there had been a greater spread of fatigue lives. Neverthe- 
less the tests appeared to be valid and do provide a measure of 
the influence of biaxial ratio on the fatigue life. 
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Results 

The results o£ the various types o£ tests and specimen 
parameter variations have been collected together in tabular 
£orm as shown in Tables 2-5. These tables illustrate the 
results o£ similar tests on each laminate and how the 
mechanical property examined was in£luenced by various 
material and test parameters. The parameters examined were 
hole diameter, biaxial load £actor, thickness, and type o£ 
specimen modi£ication. The mechanical properties examined 
were the static or ultimate strength, the £atigue li£e, and 
the residual strength. The details o£ each property varia- 
tion will be discussed on the basis o£ its variation illus- 
trated in the accompanying tables. 

Static Strength Variations 

Although it was not a primary part o£ this study, an 
initial series o£ static tests was per£ormed in which the 
hole diameter was varied £rom 12.7 to 50.8 mm to determine 
the signi£icance o£ any hole diameter in£luence on the static 
strength. These tests were per£ormed on the quasi-isotropic 
con£igurations in both the eight and sixteen ply thicknesses. 
Thus these tests were expected to provide indications o£ 
both hole diameter and specimen thickness in£luences on the 
static strength o£ these materials. 

All static and residual strength calculations, as well 
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as the applied stress level determinations, are based on 
the applied force divided by the net section area 
The net section width is taken to be the distance between 
the tips of the stress equalization slots minus the diameter 
of the center hole. The stress concentration factors associated 
with the holes were ignored since they were the same for all 
tests and would thus not affect any of the comparisons. An 
example of the stress calculations is included here so that 
there will be no misunderstanding of the implications of the 
subsequent comparisons. 

Specimen 5-1. 

Uniaxial tension: k = 0 

Hole diameter: 12.75 mm 

Net section area: 

CA^g^): 171.20 - 12.75 = 158.45 X 1.11 = 175.88 mm ^ 

Failure load (F^) : 64 ,500 N 

F 

Critical stress: = 567 MPa . 

NET 

The influence of hole diameter variation on the static 
strength of [0,±45,90]^ laminates is shown in Table 2. The 
biaxial load ratio, k, was equal to zero for all of these 
tests. It is seen that the static strength decreases 
slightly as D increases from 12.7 to 25.4 mm and then in- 
creases as D is increased to 50.8 mm. However, the entire 
variation is less than 20 percent so the effect does not 
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appear to be o£ extreme importance. It is also noted that 
the two identical tests 5-3 and 23-1, which were conducted 
from different specimen lots and over a time interval of 
many months, provided results that were nearly identical. 

A similar series of tests was performed on the 16 ply 
quasi- isotropic laminate, [0,±45,90]„ and the results are 
seen in Table 3. The biaxial load ratio, k was equal to 1.0 
for all of these tests. The same pattern seen in Table 2 
also exists in Table 3. The static strength decreased slightly 
and then increased, although the breaking strengths for D = 12.7 
and D = 50.8 mm were approximately equal in this case. If the 
values for the identical tests 9-2 and 9-3 are averaged, the 
maximum variation is less than ten percent. The average 
strength values from Tables 2 and 3 are 375 and 454 MPa, 
respectively. While this represents a significant difference, 
it cannot be attributed to any single effect because both the 
thickness and the biaxiality are different. It is suspected 
that most of the difference is due to the biaxial effect but 
there is not adequate data to verify this conclusion. 

Since the principal thrust of this research program was 
the determination of biaxial load effects, two test series 
provided data about the biaxial effect as seen in Table 4. 

These tests exhibit a significant influence of applied load 
biaxiality, for the [0,±45,90]^ laminates, while the results 
are negligible for the [0,±45,0]2g tests. For the eight ply 
quasi-isotropic laminates. Table 4 shows that the residual 
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strength increases approximately 50 percent as k increases 
from 0 to 1,0. These results coincided with a change in 
fracture path from perpendicular to the primary load axis 
to approximately 45° to it. However, it is not suggested 
that the resolved area be used to reduce the magnitude of ; 
the biaxial effect because it has not yet been possible to 
examine this matter in adequate detail. 

A pair of tests which provided an indication of the 

biaxial effect for the [0,±45,0]„ laminates are also in- 

/ s 

eluded in Table 4, The k variation was limited to the 
range 0 to 0.5 because of the weakness of this laminate in 
the 90° orientation. It is evident that no biaxial effect 
exists in these two tests, but additional tests are needed 
to confirm or disprove this behavior. 

Fatigue Tests 

A number of fatigue tests have been performed to 
ascertain the possible extent of cyclic load biaxiality on 
the fatigue life and residual strength. For these tests the 
specimens were subjected to in-phase cyclic loads at a fre- 
quench of lOHz until they failed, or until 10^ cycles were 
sustained. Any specimens lasting for 10^ cycles were then 
tested for their residual strength. The results of both the 
fatigue life and residual strength tests are seen in Table 5. 
In general the objective was to apply cyclic loads at levels 
of 70-80 percent of the static ultimate strength. However, 
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the results given in Table 5 indicate that an extremely 
broad range of cyclic lives were found to exist. 

The results of the test series on the [0,±45,0]2g 
laminates indicate a very strong and detrimental influence 
of load biaxiality on the fatigue life. However, as noted 
in the table. Specimens 11-5 and 11-4 did not fail through 
the center hole but failed due to the breaking of the loading 
tabs. Thus it seems evident that the biaxial effect would 
not have been nearly as great if the premature tab failures 
had not occurred. Because of this difficulty with tab 
failures, the series of four specimens, 7-1, 7-3, 7-4 and 
7-5, was devoted to specimen modifications intended to enhance 
specimen failures through the center hole. These attempts 
were only marginally successful so they are not discussed in 
this section. 

The results of a series of fatigue tests on [0,145,90]^ 
laminates, in which only the biaxial load ratio is varied, 
are also presented in Table 5. These tests showed that both 
the uniaxial and equal biaxial tests resulted in relatively 
short fatigue lives. However, the test with k = 0.5 resulted 
in a cyclic life approximately ten times the k = 0 and k = 1.0 
results. Because of this behavior an additional test, 23-5, 
was performed with k = 0.5, which displayed even greater 
resistance to fatigue failure. The appearance of the broken 
specimens suggested that these results may represent a valid 
biaxial effect. The k = 0 and k = 1.0 specimens both failed 



25 


in a rather brittle fashion along planes parallel to one of 
the layers (the 90° and 45° directions). However, the 
specimens with k = 0.5 failed at an angle between the 90° 
and 45° directions only after a very extensive amount of 
damage had been accumulated. It was thus concluded that 
when the preferred failure plane coincided with a ply 
orientation the failure could occur with relative ease. 
However, when failure along another plane would be anticipated 
a considerable amount of damage accumulation is necessary, 
thus resulting in a longer fatigue life. These tentative 
conclusions are also in need of further verification. 

The results of the last series of fatigue tests on 
[0,±45,0] laminates are also shown in Table 5. As discussed 
previously, all of these specimens had been modified as 
shown in Fig. 5. This modification was successful in 
eliminating the problem of tab failure, but the results were 
still not without inconsistency. 

It is seen from Table 5 that the general trend of these 
fatigue lives was toward increasing fatigue lives with in- 
creasing biaxiality in the range 0 < k < 0.25. However, the 
differences between tests 24-2 and 24-3 on one hand and 
24-4, 24-5 and 24-6 are immense, being more than three orders 
of magnitude in the fatigue life without taking into consider- 
ation the considerable residual strength. It is possible 
that the residual strength of test 24-3 was improperly re- 
corded, since it was obtained only from a peak-reading 
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digital indicator which has sometimes given erratic readings 
in the past. The conclusion which may be validly drawn from 
these tests is that the fatigue life does, in fact, in- 
crease with increasing biaxiality in the range tested. The 
fractured specimens displayed the same trend as the 23 
series, in which the amount of damage accumulated by the 
specimens was much greater when the failure could not occur 
along one of the fiber directions. Thus, it is expected 
that the fatigue life would begin to decrease as the biaxial 
factors increased beyond 0.25. 
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CONCLUSIONS 

A combined theoretical and experimental research program 
has been performed to examine the influence of biaxial applied 
loads on the mechanical properties of composite materials. 

The theoretical program consisted of developing a multi- 
layer finite-element program for examining the stress and 
damage distributions in composite laminates subjected to 
biaxial applied loads. The program was developed and 
operated with the following results. 

(1) The initial stress distribution normal to the fiber 
direction, but in the plane of the layer, near the 
hole determines the location of the initial damage 
zone and the applied stress necessary to cause 
failure. 

(2) The initial stress distribution in the fiber direc- 
tion near the hole predicts the failure path. 

(3) No delamination was predicted although it does occur 
early in the actual failure process. 

Five series of biaxial fracture toughness tests, comprised 
of 30 specimens, have been performed with the following con- 
clusions . 

(1) The size of the center hole did not have a strong 
influence on the failure load. 
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(2) Applied load biaxiality had a signficant influence 
on the static strength of [0,±45,90]^g laminates 
but a negligible influence on the [0,±45,0]^g 
laminates . 

(3) Biaxial cyclic loads had a very significant in- 
fluence on the cyclic life, although in some cases 
the influence was positive while others it was in- 
consistent or negative. The difference in behavior 
was determined by the manner in which the biaxial 
load altered the direction of the fracture path. 
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TABLE 1. Characteristics and Test Conditions o£ 
5208-T300 Graphite/Epoxy Biaxial Specimens 



TEST 

SERIES 

SPECIMEN 

NUMBER 

LAMINA 

ORIENTATIONS 

BIAXIALITY 

RATIO 

k 

HOLE 

DIAMETER 
D (mm) 

TYPE OF 
TEST* 


5-1 

[0,±45,90]^ 

0 

12.75 

1 

5-2 

ti s 

If 

50.98 

tt 

9-4 

[0, ±45,90)23 

1.0 

tt 

It 

9-5 

If 

12.75 

tf 

9-2 

[0, ±45, 90)23 

1.0 

25.53 

1 

9-3 

f f 

tt 

It 

5-3 

[0,±45,90) 

0 

tt 

tt 

5-4 

M S 

0.50 

tt 

tt 

5-5 

tt 

1.0 

tt 

tt 

11-1 

[0, ±45,0123 

0 

25.55 

tt 

11-2 

0.50 

tt 

tt 

11-3 

[0, ±45, 0123 

0 

25.48 

3 

11-4 

0.50 

25.53 

2 

11-5 

II 

0.25 

tf 

tt 

7-1 

[0,±45,0) 

f f 

25.55 

3 

7-3 

II s 

ff 

50.80 

2 

7-4 

M 

It 

25.53 

tt 

7-5 

II 

II 

50.93 

tt 

23-1 

[0,±45,90) 

0 

25.47 

1 

23-2 

ff s 

ff 

25.45 

2 

23-3 

II 

f f 

50.80 

tt 

23-4 

II 

0.50 

tt 

tt 

23-5 

II 

tt 

tt 

3 

23-6 

11 

1.0 

If 

2 

24-1 

[0,±45,0) 

0.25 

50.80 

4 

24-2 

. II s 

tt 

ff 

3 

24-3 

If 

tt 

tt 

tt 

24-4 

It 

ff 

It 

2 

24-5 

II 

0.10 

ff 

ff 

24.6 

II 

0 

If 

tt 


1 - Static strength test at low load rate. 

2 - Fatigue life test. 

3 - Residual strength test after 10^ fatigue cycles. 

4 - Static strength test at high load rate. 
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TABLE 2. Effect of Center Hole Diameter on the Static 
Strength of [0,±45,90] Laminates 

S 


SPECIMEN 

NUMBER 

BIAXIALITY 

RATIO 

HOLE 

DIAMETER 

(mm) 

BREAKING 
LOAD ^ 
(N X 10-^) 

STATIC 
STRENGTH 
CMP a) 

5-1 

0 

12.75 

64.5 

366 

5-3 

If 

25.53 

55.2 

349 

23-1 

II 

25.45 

56.9 

356 

5-2 

It 

50.98 

55.6 

429 


TABLE 3. Effect of Center Hole Diameter on the Static 
Strength of [0,±45,90]2g Laminates 


SPECIMEN 

NUMBER 

BIAXIALITY 

RATIO 

HOLE 

DIAMETER 

(mm) 

BREAKING 
LOAD , 
CN X 10^) 

STATIC 

STRENGTH 

(MPa) 

9-5 

1.0 

12.75 

161 

473 

9-2 

It 

25.53 

133 

405 

9-3 

It 

25.53 

144 

439 

9-4 

II 

50.98 

128 

472 
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TABLE 4. Effect of Applied Load Biaxiality on the Static 
Strength of [0,±45,90] and [0,±45,90], Laminates 

S ^ S 


SPECIMEN 

NUMBER 

BIAXIALITY 

RATIO 

HOLE 

DIAMETER 

(mm) 

BREAKING 
LOAD , 
(N X 10^) 

STATIC 

STRENGTH 

(MPa) 


[0,±45, 

90] g Laminates 


5-3 

0 

25.53 

55.2 

349 

23-1 

ff 

25.45 

56.9 

356 

5-4 

0.5 

25.53 

64.9 

405 

5-5 

1.0 

M 

81.0 

507 


[0,±45, 

0]2g Laminates 


11-1 

0 

25.55 

161 

505 

11-2 

0.5 

tl 

162 

505 
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TABLE 5. Effect of Cyclic Load Biaxiality on the Fatigue 
Life and Residual Strength of Composite Laminates 



SPECIMEN 

NUMBER 

BIAXIALITY 

RATIO 

CYCLIC 

STRESS 

LEVEL 

(MPa) 

FATIGUE 
LIFE , 

(Cycles X IQ-^) 

RESIDUAL 

STRENGTH 

(MPa) 



[0,±45,0] 

» Laminates 
L s 



11-3 

0 

111 

1000 

514 

11-5* 

0.25 

370 

442 

— 

11-4* 

0.50 

378 

23 

— 


[0,±45,90]g Laminates 



23-3 

0 

338 

29 

— 

23-4 

0.5 

339 

447 

— 

23-5* 

11 

341 

1000 

416 

23-6 

1.0 

343 

42 

— 


[0,±45,0] 

Laminates 

S 



24-6 

0 

335 

0.15 

— 

24-5 

0.10 

338 

0.50 

— 

24-4 

0.25 

356 

0.44 

— 

24-3 

11 

332 

1000 

451 

24-2 

ft 

305 

1000 

392 


Specimen did not fail through center hole 
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Fig. 1. Stress components on the L,T,z coordinate system. 
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Figure 2. Biaxial test system. 
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Figure 3. Biaxial specimen design. 
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Fig 


Photograph o£ typical specimen prior to testing. 



Fig. 5. Modification of specimen to reduce tab failures 
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(b) 

Fig. 6. Close-up photograph of specimen 9-2 showing, (a) 
prior damage and (b) final fracture condition. 
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Fig. 7. Photograpli of fractured specimen 7-5 showing other 
modifications . 
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ABSTRACT 


A three dimensional finite-element computer program 
has been developed to analyze layered fiber-reinforced 
composite laminate. This program is capable of: (1) cal- 

culating the detailed stress distribution, (2) identifying 
the damage zone and mode of failure, (3) analyzing the 
damage accumulation, and (4) determining the ultimate 
strength of the composite laminate. 
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1. INTRODUCTION 


The high strength and light weight properties of 
composite materials are causing them to assume increasing 
structural applications in all weight-sensitive structures, 
especially aircraft and aerospace structures. These same 
properties are also leading to increased utilization in 
other areas of the transportation sector, such as automobiles, 
where the reduction of weight leads to greater gasoline 
economy. Thus, as the structural applications of composite 
materials increase, the necessity of knowing more about their 
mechanical properties, especially their failure character- 
istics becomes more urgent. 

In the past two decades, together with the creations of 
new composite materials and the developments of new fabrica- 
tion techniques, various aspects of physical properties of 
composite materials have been investigated and different 
approaches for the better understanding of the material 
behaviors have been proposed. These efforts are published 
in the Journal of Composite Materials, ASTM proceedings, 
textbooks, and other periodicals. It is worthwhile to 
mention the work of Wendt et al. [1], Calcote [2], Garg et 
al. [3], Chamis [4], and Jones [5]. 

Since composite materials are heterogeneous and 
anisotropic, the analysis of small-scale specimens or 
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full-size structural elements made of composites is much 
more difficult than that of homogeneous and isotropic 
materials. Only in the most recent years, the development 
of techniques in numerical analysis together with large 
capacity and high speed computers make it possible to 
analyze layered fiber-reinforced composite laminates in 
detail [6-9] . 

Lee performed a three-dimensional stress analysis for 
a biaxially loaded composite laminate (with a centered 
hole) , which consists of several fiber-reinforced composite 
layers each with a specified fiber orientation, by finite 
element method [10]. In this work, a more sophisticated 
computer program was developed to analyze the damage 
accumulation including progressive failure of a layered 
fiber-reinforced composite laminate and to determine the 
ultimate strength of the composite laminate. This computer 
program is written in STAR-FORTRAN language and was operated 
on the CDC STAR-100 computer system at the NASA-Langley 
Research Center, Hampton, Virginia. 
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2. PROBLEM DESCRIPTION 

Let the specimen under consideration be a rectangular 
plate made o£ layered composite material, with a centered 
circular hole, subjected to in-plane biaxial loading (Fig. 

1) . Each layer is a unidirectional fiber reinforced com- 
posite and the direction of the fiber is characterized by 
the angle ^ between the fiber and x-axis. Thus, the 
arrangement of n layers for the entire laminate can be 
represented by ><J> 2 > • • • • • In this paper, attention 

is focused on a symmetric ply configuration which may be 
denoted as: 

[<f>l»<f>2>‘‘‘»‘l’i]s “ ' (2.1) 

The top four plies of a symmetric ply configuration 
[0,45,-45,90]^ is shown in Fig. 2. 

It is noted that the symmetric configuration has a 
mirror symmetry with respect to the x-y plane at z=0 and 
a periodicity of 180° about z-axis [10]. If the loading 
condition also has the same characteristics, only a quarter 
of the entire laminate needs to be analyzed. 
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Two basic assumptions are made in this analysis: 

1. Each layer is modeled as a linear elastic continuum 
with transverse anisotropy about fiber orientation; 
individual fibers and matrix are not considered; 

2. Unless delamination occurs, the layers are assumed 
to be perfectly bonded together; no adhesive layer 
is included in this study. 
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3. STRESS-STRAIN- RELATION 


The general stress-strain relation for a linear elastic 
anisotropic solid can be written as [11] : 



S • ■ 0 • 

ijmn mn ’ 


(3.1) 


where 


S. . 

ijmn 


= S. . 

J imn 


= S, 


ijnm 


= S_ 


mnij 


(3.2) 


The number of independent elastic constants is 21. Material 
symmetry, when it exists, imposes restrictions on the 
material constants, thus further reducing the number of in- 
dependent elastic constants. Since, for each layer, the 
unidirectional fiber reinforced composite has a axis of 
symmetry about the fiber orientation, the number of inde- 
pendent elastic constants is reduced to five and, in engi- 
neering notation, the stress-strain relation could be 
expressed as: 
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E~^ 

c-1 

■^LT^L 

c-1 

■^LT^L 

0 

0 

0 



C fjt 



e” ^ 

r-l 

■^Tz% 

0 

0 

0 


a? 


= 


■^Tz^T^ 

e'^ 

0 

0 

0 


^z 

^Tz 


0 

0 

0 

*^Tz 

0 

0 


^Tz 

^zL 


0 

0 

0 

0 

^LT 

0 


*^zL 

^LT 


0 

0 

0 

0 

0 

^LT 


^LT 

where = 0 . 5 E 

1 ^ 

■p/ (1 + '^'Pz 

) and L- 

axis 

denotes 

the fiber 


orientation. The (L,T,z) coordinate system is shown in 
Fig. 3. Now it is straightforward to obtain 




<^11 

^12 

^^13 

0 

0 

0 



0 rp 


<*12 

"^22 

^^23 

0 

0 

0 



a 

z 

= 

^13 

^23 

^3 

0 

0 

0 


^z 

^Tz 


0 

0 

0 

*^44 

0 

0 


^Tz 

^zL 


0 

0 

0 

0 

‘iss 

0 


^zL 

^LT 


0 

0 

0 

0 

0 

^66 


^LT 


(3.3) 


(3.4) 


where 



53 


*^11 ~ ’ 

^22 ^ ^33 ^ ’ 



^12 

= ^13 

= E^(l 

. VtP v^^/a 

9 




^^23 

= E^(v 

Tz " ^ 




( 


^44 

" %z 







^55 

= ‘iae 

^LT 

> 





A 

= 1 - 

^^Tz 

2v = t( 1 ♦ v^j.) 

Et/El 

• 


For 

a composite 

layer 

, whose fiber 

orientation is 

denoted 

by 4> , the general stress-strain relation 

is obtained as 

[10] 

; 









<1 

> 

<^12 

<3 » 

0 

“l6 



"y 



^22 


0 

^'Z6 




= 

^13 

A 

^^23 

‘'ss ** 

0 



"z 



0 

0 

0 d„ 

<5 

0 


''^yz 

^zx 


0 

0 

0 <5 


0 


"^zx 

V 


->l'6 


d36 » 

0 

^a'a 


Yxy 


(3.5) 


, (3.6) 


where 
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dll = diiCOs‘‘({) + d22sin“4) + 2(di2 + 2d^g) sin^(|)COS , 

^22 "" diiSin“(j) + d22COs“(}) + 2(di2 + 2d^^) sin^(|)COS^(|) , 

^33 " ^33 ’ 

A 

" d^^cos^cj) + d^^sin^c}) , 

d^^ = d^^sin^cf) + dg^cos^cf) , 

^^66 = ^66 ■" <^^11 ^22 ■ 2^12 ' ^^^66^ sin^cljcos^ct, , 

A 

di2 = di2 + (^11 ’*' ^22 ' ^*^12 ' "^^66^ sin^4>cos^(f> , (3.7) 

A 

di2 = dijCOS^(f) + d22Sin^4> , 

*^23 ” dijSin^cj) + d25COS^<}) , 

dig = sin(j)COS({> [diicos^cj) - d 22 sin^<() - cos2(t)(di2 + 2dgg)] , 

d 2 g = sinc|>cos(j) [diiSin.^4) - cos 2 (j)(di 2 + 2dgg)] , 

djg = sin(f)COS({)(di2 - d23) , 

>• 

d^g = sin(j)COS({)(d33 - d^^) . 
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4. BOUNDARY CONDITIONS 

In this work, attention is focused on a special kind 
of loading condition which is indicated in Fig. 1. That 
special kind of loading condition does have a mirror symmetry 
with respect to x-y plane at z = 0 and a periodicity of 180° 
about z-axis. Thus, it is obvious to have the following 


relations 

Uj^(-x,-y,z) = -u^(x,y,z) , 
Uy(-x,-y,z) = -Uy(x,y,z) , 
u^(-x,-y,z) = u^(x,y,z) . 

Based on eqns. (4.1) and (3.6), 
derive [10] : 

(-x,-y,z) = Oy (x,y,z) , 
(-x,-y,z) = (x,y,z) , 

= a^y(x,y,z) , 

= -a^^(x,y,z) 

o-zy(-x, -y,z) = -a^y(x,y,z) 


(4.1) 


t is straightforward to 


(4.2) 


Therefore, it is only needed to analyze a quarter of the 
entire laminate, which is shown in Fig. 4. 
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Referring to Fig. 4, for any point, pt. I(x,y,z), on 

the surface Sj, there is a conjugate point, pt. J(-x,-y,z), 

on the surface Sj*, both surfaces Sj and Sj have the same 

unit normal vector (n ,n ,0). The stress vector, t, at a 

X y ^ 

point on either Sj or Sj is related to the stress tensor at 
that point and the unit normal vector of that surface as 
follows : 


ti = o. .n. , 


or, specifically. 


(4.3) 


t 

t 

t 


X 

y 

z 


°x^x ^ ^xy^y > 
V^y ^ ^xy^x » 


a n + o n . 
zx X zy y 


(4.4) 


Then the boundary conditions on surfaces Sj and Sj may be 
specified as: 


Uy(l) - -Uy(J) , ty(l) « , 

Uj(I) ■ Uj(J) , tj(I) = . 


(4.5) 
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The rest of the boundary conditions simply are 


1. On top surface, z = h , 


0 =0 =0 = 0 ; 
z zx zy ’ 


(4.6) 


2. On symmetric surface, z = 0 , 


u = 0 =0 = 0 : 

z zx zy ’ 


(4.7) 


3. On surfaces , and S 2 , x = ±W 


0 "=k 0,0 =0 =0: 

X ’ xz xy ’ 


(4.8) 


4. On surface , y = L , 


0 = 0,0 =0 =0; 

y » "^yx yz ’ 


(4.9) 


5. On surface of the centered hole, r = rQ , 


rr 


= 0 


re 


= 0 


rz 


= 0 


(4.10) 


where 0 and k will be called the applied stress and the 
biaxial load factor. The stress boundary conditions speci- 
fied in eqns. (4,5-4.10) were converted into equivalent 
nodal force boundary conditions for use in the computer 


program. 
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5. FINITE ELEMENT ANALYSIS 

In this work, the three dimensional finite element mesh 

consists of a family of 8 -node elements. A typical element 

with the nodal points being numbered for the purpose of dis- 

T 

cussion is shown in Fig. 5. Let [6] = [u^^ ,U2 , • • . >Ug] be 
the nodal point displacement vector and 

U,(i) 

[u. ] E u^(i) ) i = (5.1) 

Uji) 

The shape functions are 

N^ = (1 + C)(l - n)(l + 0/8 , 

N 2 = (1 + 5)(l + n)(l + 0/8 , 

N3 = (1 - K ) (1 - n) Cl + ^)/8 , 

N^ = (1 - ?)(l + n)Cl + 0/8 , 

(5.2) 

Ng = (1 + c) (1 - n) Cl - 0/8 , 

Ng = (1 + ?)C1 + n)Cl - 0/8 , 

N^ = Cl - 5)(i - n)Cl - 0/8 , 

Ng = (1 - K ) (1 + n) Cl - ^)/8 . 
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Then the displacement field within this element can be 
written as: 



where, for i = 1,2,. ..,8, 


''l.X 

0 

0 

0 

N. 

0 

0 

0 

a- / Az 

0 

a^/Az 

N. 

a^/Az 

0 

N. 

1 ,x 

''’i.y 

N. 

1,X 

0 


a^ = . 


(5.6) 


(5.7) 
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And 

the coordinates (x,y) 

and (5 


''i.s 


3x/35 

3y/9^ 


"i,n 


3x/9n 

3y/3n 



N. 


N. 


1,X 


1,X 



= [J] 



N. 


N. 


i,y 


i,y 


, C5.8) 


where, referring to Fig. 5, 

|(i-n) (Xj^-x2) + (i+ri) CX2-X4) (1-n) (yi'y3) + (i+n) 


[J] = 


1 

?1 


(1+^) (x2-Xj^) + (i-c) CX4-X3) ( 1 + 5 ) (y2-yi) + (i-0 (y4-y3) 


Then the element equation which links the nodal point dis- 
placement vector [6] and the nodal point force vector [F] 
can be written as [12] : 

[K] [6] = [F] , (5.9) 

where 

[K] = I [B]'^[d'][B] dV , 

V 

1 

[B]'^[d ] [B] Az det[J] d5 dq dc , (5.10) 

-1 

and [d ] is specified in eqn. (3.6). The final governing 
equation is simply the sum of all the element equations. 
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6. THE FAILURE CRITERION 

After the governing equation of the laminate has been 
solved, the components of displacement for each nodal point 
can be obtained and then the six components of stress 
[o] ^ voint within a certain 

element can be found as 

[a] = [d^][B][6] . (6.1) 

The stress components in (L,T,z) coordinate system can be 
written as 


li 

a cos^tj) + 

OySin^cj) + 

axySin2(j) , 

(6. 

,2) 

Q 

H 

II 

a sin^i}) + 

X. 

ayCos^(J) - 


(6. 

,3) 

Q 

N 

II 

“z ■ 



(6. 

.4) 

^LT 

■ (“y ■ “x) 

sin(|)COS({) 

+ o^yCOsZcf) , 

(6. 

.5) 

^Lz 

= o^^cos* ♦ 

Oy^sincj) ; 

f 

(6. 

.6) 

'^Tz 


+ cfy^coscj) 

• 

(6. 

.7) 


To identify the damage zone, attention is focused on three 
types of damage: breakage of fibers, failure of matrix, 

and separation of layers (delamination) . The failure 
criterion for these three modes is discussed as follows: 
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1. Fiber Breakage 

For this failure mode, the stresses at the center 
of each "element are taken as the representative of that 
element. An element is identified as the damage zone 
of broken fibers if either 

, ( 6 . 8 ) 


or 


^LS 



+ 0 


Lz 



'FS 


(6.9) 


In either case, the stiffness matrix of that element, 
referring to eqn. (3.4), will be reduced to zero. 

2. Matrix Failure 

For this failure mode, the stresses at the center 
of each element are taken as the representative of 
that element. An element is identified as the damage 
zone of matrix failure if either 

’ ( 6 . 10 ) 


or 


'TS 


= (o 


TL 


+ 0 


Tz 


) ^ > 0 


MS 


( 6 . 11 ) 
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In either case, some components in the stiffness matrix 
of that element are reduced to zero, i.e.. 


[d] 


*11 

0 

^13 

0 

0 

0 


0 

0 

0 

0 


*13 

0 

^33 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 

^55 

0 


0 

0 

0 

0 


( 6 . 12 ) 


which means Orj,, and of that element will be 

reduced to zero irrespective of the strain field. 


3. Delamination 


For this failure mode, the stresses at the center 
of the interface between two layers are taken as the 
representative for the interface of the two adjacent 
elements. An interface of those two adjacent elements 
is identified as the damage zone of delamination if 
either 


^z ^ 


'DN 


(6.13) 


or 


°2S = (■’Lz' * •’tz')'' 5 '’dS • 


(6.14) 
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In either case, some components in the stiffness matrix 
of those two elements are reduced to zero, i.e.. 


[d] = 


11 

^2 

0 

0 

0 

0 


12 

42 

0 

0 

0 

0 


0 0 


0 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 


( 6 . 15 ) 


which means o^, and of those two elements 

will be reduced to zero irrespective of the strain 
field. 

The failure strength, Opg, *^DN’ 

and Ojjg, should be obtained experimentally and taken 
as the input data for the computer program. Checking 
the stresses at the center of each element and the 
stresses at the center of interface of each element 
permits the identification of the damage zone and its 
mode of failure for a specified amount of applied 


stress . 
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7 . THE PROCEDURES 

To begin with, fix the biaxial load factor k and set 
the applied stress a to be unit. After the governing 
equation of the laminate has been solved, it is straight- 
forward to find the representative stresses for each element 
and for each interface of two elements of adjacent layers, ■ 

af, <^LS’ *^T’ ^TS’ ^zS' Define six ratios R^,R 2 »...,R^ 

as follows: 

Rl = , 

R 4 = ’ 

^6 = ^zS^^DS • 

Then, for each element and its associated interface, the 
damage ratio, R(Ij , is defined as 

R(I) = Max . (Rj^ , R 2 , Rj > R 4 > Rg > Rg ) > (7.2) 

where I is the number of the element. The maximum damage 

ratio, R„„^, for the entire laminate is obtained as 
in 3 .x 

\ax ^ Max.(R(I); I = 1 , 2 , 3 , . . . , NE) , (7.3) 
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where NE is the number o£ elements. Therefore, by setting 
the applied stress a = 1/R„„^, the computer program can 
tell the location of the initial damage zone and its mode 
of failure. For the analysis of the damage accumulation, 
the subsequent procedures could be described as follows: 

Step 1: For a given set of input data, namely, the 

applied stress a, the biaxial load factor k, the 
locations and the modes of failure of the damage zones, 
the element equation for each element, eqns . (5.9,5.10), 
and the final governing equation for the entire laminate 
will be obtained. Of course, the local stiffness matrix 
d^j of the damage zones should have been modified 
according to the modes of failure as indicated in the 
previous section. 

Step 2: Solve the governing equation for the entire 

laminate and obtain the stress distribution for each 
element and its associated interface. For each element 
and its associated interface, compare the representative 
stresses with the failure strength and see whether any 
damage happens; if yes, change the stiffness matrix 
accordingly; otherwise, calculate the damage ratio, 

R(I). 
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Step 3: Going through all the elements and their 

associated interfaces, if the damage zones are enlarged, 
keep the applied stress unchanged together with the 
enlarged damage zones and modes of failure as input 
data and go back to Step 1. Otherv>rise, calculate the 
maximum damage ratio, and increase the applied 

stress by the factor SO back to Step 1. 

It is seen that there are two possibilities for a given 
applied stress, damage zones, and modes of failure The 
first possibility is that the damage zones will not be 
enlarged and then that kind of damage is said to be localized 
at the applied stress specified. The second possibility is 
that the damage zones will be enlarged at a constant applied 
stress until the entire laminate is broken into pieces. In 
the second case, it is said that progressive failure occurs 
and the ultimate strength of the composite laminate is 
identified. 
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8. NUMERICAL RESULTS 

First, let the input data relevant to this work be 
listed as follows: 

1. Material constants: 

Ej^ = 134.58 X 10^ MPa , 

E.J. = 9.68 X 10^ MPa , 

= 0.264 , (8.1) 

^Tz 0«435 , 

= 4.48 X 10^ MPa . 

2. Failure strength: 

Op^ = 1241 MPa , 

Opg = 82.7 MPa , 

*^MN ~ 55.2 MP a , 

( 8 . 2 ) 

^MS ~ 82.7 MPa , 

Odn " 55.2 MPa , 

^DS ~ 89.6 MPa . 

3. Ply configuration: 

[0,45,-45,90]g = 


[90,-45,45,0,0,45,-45,90] 


(8.3) 
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4. Geometrical parameters: 

2L = 152.4 mm , 

2W = 152.4 mm , 

(8.4) 

2h = 1.07 mm , 
rQ = 12.7 mm . 

5. Finite element mesh: 

No. of nodal points =1750 , 

No. of degrees of freedom = 5250 , 

(8.5) 

No. of elements = 1248 , 

3 

Volume of smallest element = 0.175 mm 

For illustrative purpose, the numerical results of 
three cases, k = 0,0. 5,1, are presented as follows: 

Case 1: uniaxial loading, k = 0 

The maximum tensile stress along the fiber among 
the top four layers, max. is plotted in Fig. 6 as 

a function of 0 at r = 13.1 mm. It is noticed that 
the overall max. Oj^ occurs at 0 = 0° (180®) in the 
fourth layer in which the fibers are parallel to the 
y-axis--the loading direction and the maximum stress 
concentration factor is as high as 7.92, which means 
the fourth layer carries the major portion of the load. 
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The maximum normal stress perpendicular to the fiber 
among the top four layers, max. o^, is plotted in 
Fig. 7 as a function of 0 at r = 13.1 mm. It is 
noticed that the overall max. (= 0.56a) occurs at 
0 = 172.5° (- 7.5°) in the first layer, which indicates 
the first damage zone of matrix failure will be initiated 
in the element of the first layer centered at r = 13.1 mm 
and 0 = 172.5° (- 7.5°) when the applied stress a 
reaches 98.6 MPa (Fig. 8). At a = 138 MPa, the damage 
zone is shown in Fig. 9. However, those damages are 
localized, which means that the damage zone will not be 
enlarged if the applied stress is not increased further. 
At a = 156 MPa, the first fiber breakage occurs in the 
element of the fourth layer centered at r = 13.1 mm 
and 0=0° (180°) as shown in Fig. 10 and that leads to 
progressive failure as indicated by Fig. 11 and Fig. 12. 
In this case, the ultimate strength of the composite 
laminate is identified to be 156 MPa. 

Case 2; Biaxial loading, k = 0.5 

Similarly, max. and max. a^ are plotted as 
functions of 0 at r = 13.1 mm in Fig. 13 and Fig. 14, 
respectively. Again, the overall max. Oj^ occurs at 
0=0° (180°) in the fourth layer but the maximum stress 
concentration factor is reduced to 6.57 due to the load 
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applied parallel to the x^axis. The overall max. 
o,j. (= 0.466a) occurs at 0 = 172.5° (-7.5°) in the first 
layer, which means the first matrix failure will happen 
when the applied stress o reaches 119 MPa (Fig. 15). 

The damage zone at a = 175 MPa is shown in Fig. 16. 

At a = 189 MPa, fiber breakage occurs and that leads 
to progressive failure as indicated by Figs. 17,18; 
thus, the ultimate strength in this case is obtained to 
be 189 MPa. 

Case 5: Equal biaxial loading, k = 1 

The maximum tensile stress along the fiber among 
the top four layers, max. plotted in Fig. 19 as a 
function of 6 , is noticed to have four peak values - 
5.226a, 5.252a, 5.214a, 5.245a - at 6 = 0° (180°), 

45° (225°), 90° (270°), 135° (315°), respectively. This 
means that, first, the maximum stress concentration ( 
factor is further reduced to the neighborhood of 5.2 
and, second, the final broken path of the specimen may 
be one of the four lines: 6 = 0°-180°, 45°-225°, 

90°-270°, 135°-315°. Figure 20, which shows max. a^ 
as a function of 0, indicates that the overall max. 
a^ (= 0.382a) occurs at 0 = 165° (-15°) in the first 
layer. This means the first matrix failure will happen 
when the applied stress a reaches 144 MPa (Fig. 21). 
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The damage zones at a = 172 MPa and at a = 192 MPa 
are shown in Fig. 22 and Fig. 23, respectively and 
it is noticed that those damages are localized. At 
a = 193 MPa, fiber breakage occurs and that leads to 
progressive failure as indicated in Figs. 24,25; thus, 
the ultimate strength in the equal biaxial loading 
case is obtained as 193 MPa. 
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9. DISCUSSIONS 


This three dimensional finite element computer program 
can be used to calculate the six components of stress 

^^yz’ ""zx’ V equivalently, o^, a^, o^, o^g, 
Ots> point in the composite laminate. Based on 

this stress distribution and the given failure criterion, 
the damage zones and the modes of failure can be identified. 
According to the modes of failure, the stiffness matrix of 
those elements in the damage zones will be modified and 
then subsequently the new stress distribution is obtained; 
This new stress distribution determines the enlargement of 
the damage zones or the higher load which the composite 
laminate is able to carry. This process can go on and on 
until progressive failure occurs and therefore the ultimate 
strength of the composite laminate is obtained. 

From the three cases presented in the previous section, 
it is noticed that: (1) the initial stress distribution of 

Orp near the edge of the hole determines the location of 
initial damage zone and its mode of failure as well as the 
applied stress for initial damage; (2) the initial stress 
distribution of near the edge of the hole predicts the 
final broken path of the specimen. 


It has also been noticed that, during the whole process 
of damage enlargement and progressive failure, the computer 
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program has never detected any delamination . Many authors, 
including Spilker and Chou [13], Altus et al . [14], Raju and 

Crews [15] , have investigated the free edge effects in 
composite laminate and reported that both normal stress and 
shear stress between two composite layers have singularities 
near the free edge and therefore, logically, delamination 
should happen prior to the occurrence of other modes of 
failure. This program's failure to detect the occurrence 
of delamination is most likely due to the coarseness, 
relatively speaking, of the finite element mesh near the 
edge of the hole. However, further refinement of the finite 
element mesh to the required level at this stage is practically 
impossible and, hence, additional effort and a new approach 
will be necessary in this respect. 
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y 



Fig. 3 The Stress Components in |L, T, z) Coordinate System 
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MAX. ar^/cr 



Fig. 6. MAX. or. vs 6 at r =13.1 mm.j k = 0 



Fig. 7. MAX. <r_ vs 8 at r = 13.1mm., k = 0 
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Fig. 9. Damage zone at O - 138 MPa, k - 0. Matrix Failure ^ Fiber Breakage 
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Fig. 11. Damage zone at a - 156 MPa. k - 0. Matrix Failure ^ Fiber Breakage Hj 
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MAX. o-^/cr 



Fig. 13. MAX. vs 9 at r = 13.1 mm, k = 0.5 




MAX. <T^/<T 



Fig. 14. MAX. a,j, vs 9 at r = 13.1 mm, k = 0.5 



Fig, 15. Damage zone at O - 119 MPa, k *= 0. Matrix Failure Fiber Breakage 
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Fig. 24, Damage zone at O - 193 MPa, k - 0. Matrix Failure ^ Fiber Breakage 
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